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ABSTRACT 

There has been much recent work dedicated to exploring secondary correlations in 
the mass-metallicity relation, with significant dependence on both the star formation 
rate and Hi content being demonstrated. Previously, a paucity of molecular gas data 
(combined with sample selection bias) hampered the investigation of any such relation 
with molecular gas content. In this work, we assemble a sample of 221 galaxies from a 
variety of surveys in the redshift range 0 < z < 2, to explore the connection between 
molecular gas content and metallicity. 

We explore the effect of gas mass on the mass-metallicity relation, finding that the 
offset from the relation is negatively correlated against both molecular and total gas 
mass. We then employ a principle component analysis technique to explore secondary 
dependences in the mass-metallicity relation, finding that the secondary dependence 
with gas mass is significantly stronger than with star formation rate, and as such the 
underlying ‘Fundamental Metallicity Relation’ is between stellar mass, metallicity, and 
gas mass. In particular, the metallicity dependence on SFR is simply a byproduct of 
the dependence on the molecular gas content, via the Schmidt-Kennicutt relation. 
Finally, we note that our principle component analysis finds essentially no connection 
between gas-phase metallicity and the efficiency of star formation. 

Key words: galaxies: evolution - galaxies: formation - galaxies: abundances - galax¬ 
ies: statistics 


1 INTRODUCTION 

Much of the quest to understand the processes driving 
galaxy evolution has concentrated on understanding the so- 
called ‘scaling relations’, distinctive correlations existing be¬ 
tween large-scale physical parameters of galaxies. These scal¬ 
ing relations allow insight into the mechanisms underlying 
galaxy formation by providing simple, quantitative tests for 
models - both analytical and numerical - to match. Indeed, 
the ability (or otherwise) of a model to reproduce a variety 
of observed scaling relations has become a critical metric by 
which a model’s success is judged. 

In recent years, the well-known scaling relations be¬ 
tween stellar mass and metallicity (the ‘mass-metallicity 
relation’), and stellar mass and star formation rate (the 
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‘main sequence’ of galaxy evolution) have been extended 
into a three-dimensional relation, known as the Fundamen¬ 
tal Metallicity Relation (FMR; [Mannucci et al.||2010 Lara- 
Lopez et al. 20101. Galaxies up to at least 2 ~ 2.5 lie on 


the surface defined by the FMR, on which (a) more mas¬ 
sive galaxies have higher metallicities, and (b) at a given 
stellar mass, galaxies with higher star formation rates have 
systematically lower metallicities. 

The existence of the FMR has provided an excellent 
target for models to reproduce, with a variety of models 
successfully reproducing the general form of the relation, 
including analytical models (Dave et al.||2012 Dayal et al. 

semi-analytical models | |Yates &: Kauffmann||2014p , 
and SPH/radiative transfer models l Obreja et al.|2014 |. As 
noted by |Bothwell et ah (20131, these models have primarily 
explained the correlation between stellar mass, metallicity, 
and star formation rate in terms of the global behaviour of 
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gas, which can dilute the metallicity of the ISM via pristine 
inflows, trigger metallicity-enriching bursts of star forma¬ 
tion, and remove metals from the ISM in the form of galactic 
winds. 


Based on this, Bothwell et al. (20131 conducted a study 


of 4253 local galaxies, finding that the mass-metallicity re¬ 
lation also exhibited a significant secondary dependence on 
mass of atomic hydrogen (Hi), to the extent that the ‘Hl- 
FMR’ was potentially more fundamental than the correla¬ 
tion with SFR (in terms of the scatter around the respective 
relations). As with the SFR-FMR before, this has proved 
a fertile area for models, with several authors concluding 
that it is indeed the gas content of galaxies that funda¬ 
mentally underlies and drives the mass-metallicity relation. 
Dave et al. (20121 present an analytical model of galaxy 


evolution, which predicts that metallicity will anti-correlate 
with Hi mass. [Zahid et al.||20i^ present a model whereby 
the metallicity of galaxies is regulated by the gas-to-stellar 
mass ratio, such that the mass-metallicity relation can be 
seen as a consequence of the underlying metallicity/gas-to- 
stellar mass ratio scaling. [Ascasibar et al.]|2014| conduct a 
similar investigation, extending to resolved studies of local 


galaxies, and reaching similar conclusions. (See also Magrini 


et al.|20T2 Lara-Lopez k, Hopkins|2014[ Ho et al.|2015 1 for 

additional work). 

One potential surprise, however, was the lack of a clear 


dependence with molecular gas. While Bothwell et al. (20131 


found a close and tight relationship between stellar mass, 
metallicity, and Hi mass, any such relation with H 2 mass 


remained elusive. As discussed by Bothwell et al. (20131, 


this was potentially due to two factors: (1) sample selec¬ 
tion bias, resulting in a small number of galaxies covering 
insufficient parameter space, and (2) the fact that measur¬ 
ing H 2 mass requires a ‘conversion factor’ (to convert from 
its observable tracer, ^^CO), which is itself a function of 
metallicity. This lack of connection between molecular gas 
mass and the mass-metallicity relation was particularly sur¬ 
prising, because the physical processes thought to be driv¬ 
ing the Hl-FMR (including star formation-triggering inflows, 
and metal rich outflows) are intimately connected with the 
molecular gas component of the ISM. Furthermore, it has 
long been established that molecular gas and star formation 
are tightly correlated across a wide range of galaxy types, via 


the Kennicutt-Schmidt relation (Schmidt 1959 Kennicutt 


1998); given this relation, a H 2 -FMR should be expected on 


the basis of the existence of the SFR-FMR alone. As such. 


Bothwell et al. (2013) concluded that while a H 2 -FMR was 


likely to exist, quantifying its existence was beyond the scope 
of their data. 

In this work, we assemble a large sample of galaxies, 
from dwarf galaxies in the local Universe to extreme star- 
bursts at 2 ; ~ 2, to explore the connection between the 
molecular gas content of galaxies, and their gas-phase metal¬ 
licity. First, we examine the effect of molecular gas mass on 
metallicity scaling relations, by searching for a ‘H 2 -FMR’ ef¬ 
fect, whereby galaxies show offsets from the mass-metallicity 
relation which correlate with molecular (and total) gas mass. 
We then adopt a non-parametric statistical approach, and 
use a principle component analysis technique to examine 
correlations in stellar mass/metallicity/gas mass parame¬ 
ter space, demonstrating that (a) the most signihcant ‘sec¬ 
ondary correlation’ in the mass-metallicity relation is indeed 


with gas content, and (b) there is very little correlation be¬ 
tween metallicity and star formation efficiency. 

Throughout this work we use a Kroupa (2001) IMF, 
and a cosmology following the latest Planck results, with 


Ho = (67.8±0.9) km/s/Mpc and Qm ~ 0.308T0.012 (Planck 
Collaboration et al.|201^. 


2 SAMPLE SELECTION 


In this work we will be characterising the scaling rela¬ 
tions between a variety of parameters: stellar mass, metal¬ 
licity, star formation rate, and gas content (both atomic and 
molecular). We have selected our samples based on the si¬ 
multaneous availability of all of these parameters. Stellar 
masses, star formation rates, atomic gas masses, and metal- 
licities (typically measured via optical strong line diagnos¬ 
tics) are comparably easily measured, and are widely avail¬ 
able in the literature. Molecular gas masses (measured using 
observations of ^^CO) are less widely available, and it is the 
availability of this latter parameter that typically represents 
the limiting factor in our sample selection. 

In the local Universe, we use three recent surveys for 
molecular gas which fulfil these criteria excellently: COLD 
GASS (|Saintonge et al.|2011[ ), the Herschel Reference Survey 
( Boselli et al.|2014 ), and ALLSMOG (Bothwell et al. 2014a). 
Both COLD GASS and ALLSMOG are selected from the 
SDSS spectroscopic survey, and as such have available op¬ 
tical spectra for deriving metallicities. Optical spectroscopy 
for the Herschel Reference Survey (HRS) is available from 
[Boselli et al.| ( [20T^ . In order to gather better statistics at the 
low-mass, low-metallicity end of the distribution, we also use 
galaxies from the volume-limited Local Volume Legacy sur¬ 
vey, which (being a complete volume-limited sample within 
11 Mpc) contains galaxies of lower mass than appear in other 
samples. 

We also wish to extend our study beyond the local Uni¬ 
verse: in order to compare the samples discussed above to 
those at higher redshifts, we gather samples of both ‘normal’ 
star forming galaxies (selected using optical colour cuts, and 
lying on or around the main sequence) and luminous sub¬ 
millimetre galaxies, SMGs (selected via their sum-mm flux, 
and having extreme star formation rates significantly above 
the star-forming ‘main sequence’). These high -2 galaxies 
typically lie around 2 ~ 2. 

Below, we briefly describe these various samples. A his¬ 
togram of stellar masses for galaxies used in this work (sepa¬ 
rated by sample) is shown in Fig.[^ Our combined samples 
provide a mass coverage of approximately three orders of 
magnitude, from 8.5 < logM*/MQ < 11.5. In total, our 
combined sample consists of 221 galaxies from 0 < z ^ 2. 


2.1 Local Universe surveys 

2.1.1 ALLSMOG 

ALLSMOG, the Apex Low-redshift Legacy Survey for 
Molecular Gas, is an APEX large programme designed to 
survey molecular gas in local galaxies with stellar masses 
8.5 < \og[Mt./M q) < 10.0. The survey covers the local Uni¬ 
verse, 0.010 < 2 < 0.025 (= 40 <; D/Mpc <1 110), excluding 
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Figure 1. Histogram of the stellar masses of galaxies in the sam¬ 
ples used in this work. Only galaxies that obey the selection crite¬ 
ria detailed in [j^are shown. The main panel shows the combined 
stellar mass histogram, while the smaller inset panels show his¬ 
tograms of the stellar mass for each separate sample. 


The ALLSMOG galaxies used in this analysis have a 
metallicity range [8.52 <12-1- log(0/H) < 9.17], a SFR 
range [0.04 < SFR/Mq yr“^ < 12.0], and lie between stellar 
masses [8.5 < logM*/MQ < 10.0]. 

All ALLSMOG galaxies published in the first data re¬ 
lease are used in this analysis: a total of 42 galaxies. 


2.1.2 COLD GASS 


The COLD GASS survey is a large IRAM program, designed 
to survey molecular gas in massive a ~ 0 galaxies, with 
log(M*/MQ) > 10. The survey covers the redshift range 
0.025 < 2 < 0.05 (= 110 < D/Mpc < 220), and is selected 
at random from the parent GASS sample jCatinella et al. 
20101. The COLD GASS sample is discussed in detail by 
Saintonge et al. ]( |MT| ). 

As with the ALLSMOG sample, COLD GASS is drawn 
from the SDSS spectroscopic survey and has associated 
available spectra, which supply both stellar masses and star 
formation rates (as detailed in the ALLSMOG description 
above). We remove AGN (which contaminate the optical line 
fluxes, making the resultant metallically estimates inaccu¬ 
rate) from the COLD GASS sample using the BPT cut given 
by Kauffmann et al)] ( |2003| ). The COLD GASS galaxies used 
here have a metallicity range [8.87 < 12-|-log(0/H) < 9.29], 
a SFR range [0.04 < SFR/Mq yr“^ < 27], and lie between 
stellar masses [10.0 < logM*/MQ < 11.5]. 

After discarding AGN, and galaxies with discrepancies 


between their R23 and [Nii]/Hq metallicities (see [2.3.1 be¬ 
low), we are left with 115 COLD GASS galaxies. 


the very nearby Universe to ensure that the 27” APEX beam 
recovers the total CO flux. 

Upon completion, the ALLSMOG survey will have CO 
measurements (both detections and upper limits) for ~ 100 
low mass (< 10 ^°Mq) galaxies. In this work, we use the 
‘halfway’ data release presented by Bothwell et al. (20141 
which consists of 42 galaxies. Data observations and reduc¬ 
tion were undertaken as described in that work. Being drawn 
from the SDSS spectroscopic survey, all ALLSMOG galaxies 
have available optical spectra. One of the aims motivating 
the ALLSMOG survey was to study the connection between 
metallicity and molecular gas content, and so all ALLSMOG 
galaxies have the complete suite of optical line fluxes nec¬ 
essary for deriving metallicities (see ^2.3.1| for an explana¬ 
tion of the metallicity derivation). All ALLSMOG galaxies 
are selected from the SDSS spectroscopic survey, and have 
their stellar masses and star formation rates readily avail¬ 
able. Both stellar masses and star formation rates, corrected 
for aperture effects caused by the fibre, are provided by the 
MPA-JHU grouiQ Star formation rates are derived using 
optical line fitting, as described in Brinchmann et al. (20041. 
Stellar masses are likewise measured using the fitting tech¬ 
nique described by Kauffmann et al. (20031. Uncertainties 
on these parameters are typical for SDSS galaxies - approx¬ 
imately ~ 0.3 dex and ~ 0.15 dex for SFRs and stellar 
masses, respectively. 


^ http://www.mpa-garching.mpg.de/SDSS/DR7/ 


2.1.3 Herschel Reference Survey 


The Herschel Reference Survey (HRS) is a Herschel guaran¬ 
teed time key project, designed to study a complete, volume 
limited sample of galaxies at distances 15 < D/Mpc < 25. 
Within these limits the survey is highly complete. A full 
survey description is given by |Boselli et al.| ( |2010[ ). 

Lying closer than 25 Mpc, the HRS galaxies are too 
nearby to be effectively probed by a fibre survey such as 
SDSS (the 3” SDSS fibre would cover only a tiny fraction 
of the galaxy disc). However, ancillary parameters - includ¬ 
ing both stellar masses and optical spectra - are available 
from other sources. Stellar masses for the HRS are presented 


by Cortese et al. (20121, and were derived using i-band lu¬ 
minosities, with a typical uncertainty of ~ 0.15 dex. Star 
formation rates for the HRS galaxies can be calculated from 
a combination of UV and IR fluxes (UV fluxes presented by 
Gortese et aL|20T2| and IR fluxes, presented by jBendo et al.j 
20121. Briefly, we define the FUV attenuation via the IRX’ 


parameter, IRX= log(LiR/Lpuv,obs): 


A(FUV) = -0.028A® -h 0.392X2 1.094X + 0.546, (1) 

where X is ‘IRX’ ( [Burgarella et al.|2005p . SFRs are then 
calculated from the attenuation-corrected FUV luminosity: 


log SFR = log Lpuv.corr - 9.68, (2) 

which is a SFR prescription derived for use with the 
GALEX bands by Iglesias-Paramo et al. (20061. FUV,corr 
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Figure 2. Histogram of the star formation rates (left panel), molecular gas masses {central panel), and metallicities {right panel) of 
galaxies in the samples used in this work. As before, only galaxies that obey the selection criteria detailed in 0 are shown. 


is the FUV luminosity, corrected for the effects of UV at¬ 
tenuation. We have reduced their original UV luminosity to 
SFR conversion factor by a factor of 1.5, in order to convert 


quiring both molecular gas data and an optical metallicity 
measurement serves to eliminate much of the LVL sample: in 
total, 23 galaxies from LVL have the full complement of re¬ 


from their Salpeter (19551 IMF to the 

Kroupa (20011 IMF quired data. (One of these galaxies, however - UGC05364 - 

used to calculate SDSS parameters (see 

Salim et al.|2007 for has its metallicity estimated using observations of planetary 


details). 

Optical spectra for the HRS, presented by|Boselli et al.| 


120131, were taken by the 1.9m telescope at Observatoire de 
Haute Provence, and cover the spectral range 3600-6900A, 
allowing metallicities to be derived identically as for the 
SDSS samples described above. Again, we use a BPT cut 
to remove AGN from the sample, and do not use galax¬ 
ies with discrepant R23 and [Nil]/Ha metallicities (as de¬ 
scribed by |2.3.1| below). The HRS galaxies used here have 
a metallicity range [8.57 <12-1- log(0/H) < 9.05], a SFR 
range [0.70 < SFR/Mq yr“^ < 16.5], and lie between stellar 
masses [8.9 < logM,/MQ < 10.4]. 

After these cuts, we are left with 24 HRS galaxies. 


2 . 1.4 Local Volume Legacy survey 

In addition to the three dedicated CO surveys discussed 
above, we also include galaxies from the volume-limited Lo¬ 
cal Volume Legacy (LVL) survey. LVL is a IRAC and MIPS 
Legacy survey of a volume-limited sample of 258 galaxies, 
designed to be statistically complete within 11 Mpc. The 
LVL survey was designed to work in synergy with ancil¬ 
lary datasets, including the UV and Ha survey IIHUGS (11 
Mpc Ha and UV Galaxy Survey). The combined dataset is 
volume-limited out to 11 Mpc (excluding the ‘zone of avoid¬ 
ance, defined by the Galactic plane ]&] <; 20°), with a mag¬ 
nitude limit R < 15. 

Stellar masses for the LVL sample have been estimated 


by Bothwell et al. (20091, by combining photometry with a 


colour-dependent mass to light ratio, taken from the models 


of Bell & de Jong (20011. Bothwell et al. (20091 quote the 
estimated uncertainty on the stellar masses as ~ 0.3 dex. 
Gas data, both atomic and molecular, were compiled from 


a range of literature sources by Bothwell et al. (20091. Star 
formation rates for the LVL sample were derived by |Lee et al.j 
(20091, by combining UV and IR fluxes as per the method 


outlined above. The average uncertainty on these SFRs is 
approximately ~ 0.3 dex. 

The final required data product, optical metallicity 


nebulae - a unique method, not shared by any other galaxy 
used in this analysis. We therefore discard UGG05364 from 
the sample.) The LVL galaxies in our combined sample have 
a metallicity range [8.01 < 12 -|- log(0/H) < 9.31], a SFR 
range [0.04 < SFR/Mq yr“^ < 5.2], and lie between stellar 
masses [8.4 < logM*/MQ < 11.5]. 

The LVL sample used in this analysis consists of 22 
galaxies. 


2.2 Galaxies at high redshift 

2.2.1 Main sequence star-forming galaxies 

Our sample of high-z ‘main sequence’ star forming galax¬ 
ies is taken from the PHIBBS survey for molecular gas in 
1 < 2 : < 3 ‘main sequence’ galaxies ( [Tacconi et al.|[2013p . 
These ‘main sequence’ galaxies are selected via the ‘BX/BM’ 
optical colour-selection criteria (which approximately selects 
‘normal’ star forming galaxies at high- 2 ; details of this se¬ 
lection criteria are given in Steidel et al.||20()4 l. 

Galaxies in the PHIBBS survey lie at two distinct 
epochs: 2 ~ 1.2 and 2 ~ 2.2. The lower redshift galaxies are 
unsuitable for inclusion in this work, as their optical spectra 
(provided by the DEEP2 survey) do not cover any suitable 
metallicity-tracing line ratios. The galaxies in the higher red¬ 
shift bin, however, have been observed with VLT/SINFONI 
( Newman et al.||2012 l, and have publicly-available near-IR 
spectroscopy covering the Ha and [NH] lines, which we use 
to derive gas-phase metallicities (as described in ^2.3.1| be¬ 
low). 

Stellar masses, molecular gas data, and star formation 
rates for these galaxies have been taken from [Tacconi et al.j 
(20131. To briefly summarise, stellar masses have been de¬ 


rived from UV-to-IR population synthesis modelling, and 
star formation rates were calculated using a combination of 
UV, optical continuum, Ha, and 24/rm continuum fluxes. 


Tacconi et al. (20131 quote a typical stellar mass uncer¬ 


data, is compiled and presented by Marble et al. (20101. Re¬ 


tainty of 0.13 dex, and a typical SFR uncertainty of 0.15 
dex. Molecular gas data for this sample were obtained us¬ 
ing the IRAM Plateau de Bure Interferometer. Due to the 
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redshift of the sample, the sources here are observed in the 
^^CO(3 —2) line - for the purposes of calculating gas masses, 
we have converted down to equivalent ^^CO(l — 0) lumi¬ 
nosities by assuming a ‘Milky Way’-like CO excitation, ap¬ 
propriate for these normal star-forming galaxies. The high- 
2 ; ‘main sequence’ galaxies in our sample have a metallic¬ 
ity range [8.61 <12-1- log(0/H) < 9.05], a SFR range 
[30 < SFR/MQyr“^ < 245], and lie between stellar masses 
[9.8 < logM./Mg < 11.4]. 

Our final sample of high -2 main sequence galaxies con¬ 
sists of 9 objects. 


2.2.2 Sub-mm galaxies (SMGs) 


We also include a sample of galaxies selected at sub-mm 
wavelengths (typically Ssso ^ 5mJy), known as sub-mm 
galaxies (SMGs). SMGs are high-redshift dusty galaxies 
with extreme SFRs ( ^ 1000 MQyr“^). The source of the 
high SFRs exhibited by SMGs remains a matter of debate, 
but it is likely that the population contains both major 
mergers and extreme, ‘maximum starburst’ disks. Either 
way, it is clear that SMGs represent an extreme population 
lying far above the star-forming ‘main sequence’ at their 
redshifts. The inclusion of SMGs in our combined sample 
will extend the dynamic range of our analysis, in terms of 
both parameter space and the range of ISM environments 
sampled. 

Our sample of SMGs lie at 2 ~ 2, and are taken from 


the sample presented by Chapman et al. (20051. Metal 


licities were derived from Keck/NIRSPEC near-IR spec¬ 
troscopy ( [Swinbank et al.]|2004 |, which covers both the Ho? 
and [Nil] emission lines (again, see |2.3.1 below for details of 
the metallicity derivation). We note that SMGs are highly 
dusty galaxies, and as such optical emission lines may not 
fully trace regions of active star formation. However in the 
absence of well-calibrated far-IR diagnostics (as described 
in, e.g., Nagao et al 2011), these optical lines remain the 
only way to estimate metallicities in SMGs. We therefore 
proceed with the metallicity calibration, and simply warn 
readers that the metallicity values for our SMGs may be un¬ 
certain. Stellar masses for these SMGs were derived from full 
population synthesis modelling by |Hainline et al.] ( |2011| , and 
star formation rates were calculated using IR fluxes (derived 


from the far IR-radio correlation) by Bothwell et al. (20131. 


Uncertainties on star formation rates and stellar masses for 


SMGs are typically of the order 0.3 dex (as given by Mainline 
et al.|20TT and Bothwell et al.|2013 l. 


Molecular gas masses were taken from the CO survey 


by Bothwell et al. (20131. As with the high -2 main sequence 


galaxies above, molecular gas masses are based on mid-J 
^^CO measurements - either ^^CO(4 — 3) or ^^CO(3 — 2). 
We have adjusted these down to the equivalent ^^CO(l — 0) 
luminosities using the SMG CO excitation derived byjBoth-j 
well et al.| ([M3|. One SMG (SMMJ030227.73-f000653.3) 


has an unusually high ]Nii]/Hq ratio; as noted by Swinbank 


et al. (20041, this is likely indicative of contamination by 


AGN activity. As such, we exclude it from our analysis. 

The SMGs in our analysis have a metallicity range 
[8.74 < 12 + log(0/H) < 9.09], a SFR range [390 < 
SFR/MQyr“^ < 2890], and lie between stellar masses 
[10.3 < logM./Mg < 11.2]. 


Excluding SMMJ030227.73-1-000653.3, our final SMG 
sample consists of 9 galaxies. 


2.3 Deriving physical parameters 


2.3.1 Deriving metallicities 

For all galaxies in this work (with the exception of some 
galaxies in the LVL sample: see below), we derive gas-phase 
metallicities using optical strong-line fluxes. Where possible, 
we take optical line fluxes from SDSS (this is applicable 
for the ALLSMOG and COLD GASS samples). Following 
Mannucci et al (2010), we selected galaxies for inclusion in 
our ALLSMOG sample by requiring a signal-to-noise ratio 
of 25 in the Ha line - this ensures that the main optical lines 
are generally detected at high enough S/N to reliably derive 
a metallicity. All galaxies in the COLD GASS sample also 
fulfil this criteria. Mean line flux errors for our SDSS sources 
are as follows: Hoerr ~ 4%; [NH]err ~ 5%; H/3err ~ 25%; 
[OH]3727err ~ 48%; [OHI]5007err ~ 28%; [OHI]4959err ~ 
31%. The faintest line detected is generally [OH]3727, and 
on average sources in our samples have [OH] 3727 detected 
with a S/N ratio of ~ 2. All lines are individually corrected 
for dust extinction, which is estimated using the Ha/HP 
ratio, and a Milky Way extinction curve. 


For ease of comparison, we follow Mannucci et al. (20101 


in estimating metallicity based on two independent diag¬ 
nostics: (1) the N[ll]/Ha ratio, and (2) the R23 parame¬ 


4958, 5007)/H/3. We then use the 
Maiolino et al.| (|2008|) to calculate 


ter, = ([OH]3727 -f [OIH 
abundance calibrations of 
two separate metallicity measurements, one for each diag¬ 
nostic. We take as our final metallicity value the mean of 
these two metallicity measurements, discarding any galaxies 
for which the methods give answers discrepant by > 0.2 dex: 
such galaxies are likely to have inaccuracies in one or both 
metallicity diagnostics, and as such are unreliable. We note 
that 9 galaxies in our ALLSMOG sample lie at 2 < 0.019, 
and thus the [OH] 3727 line does not fall within the SDSS 
spectral range. For these 9 galaxies, instead of the R23 pa¬ 
rameter (which requires the detection of the [OH] 3727 line) 
we adopt the metallicity as listed in the MPA-JHU catalogue 
(which is derived following Tremonti et al. 2004), adjusted 
for our metallicity calibration. We then, as above, take the 
mean of this value and the value derived using the N[ll]/Ha 
ratio. 

One potential source of metallicity bias in our SDSS- 
observed samples is due to aperture effects. SDSS is a fibre- 
based spectroscopic survey, and the 3” fibre may only sam¬ 
ple the central few kpc of some galaxies, if they are nearby, 
and subtend large areas on the sky. Given the existence of 
metallicity gradients, this may cause abundance estimates 
to be biased high (though this effect is likely to be small, as 
gradients in star-forming galaxies tend to be shallow - ~ 0.1 
dex/Re; [Sanchez et al.|[2014| ). [ifewley et al.| ( |2005| ) estimate 
that a ‘fibre covering fraction’ of 20% is sufficient to recover 
the metallicity of a galaxy without aperture bias. If we re¬ 
move the ALLSMOG and COLD GASS galaxies which fail 
this criteria (3 ALLSMOG galaxies, and 33 COLD GASS 
galaxies), our results are unaffected. For the remainder of 
this work, we therefore proceed with the full samples as de¬ 
scribed above. 

To recap the various sources of the non-SDSS spec- 
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troscopy, HRS galaxies have long-slit integrated spec¬ 
troscopy presented by Boselli et al. (20131, high-z ‘main se¬ 
quence’ galaxies have near-IR spectra presented by |Newman| 
et al. (|2012 I, and the SMGs have near-IR spectra presented 


by Swinbank et al. (20041. For the galaxies at high-redshift 
(both main-sequence and SMGs), the full suite of optical 
strong lines is not available for each source. Instead, we es¬ 
timate their metallicities based on the N[ll]/Ha ratio alone, 
again using the abundance calibrations provided by|MaioIino| 


et al. (20081. 


As mentioned in |2.1.4| above, galaxies in the LVL sam¬ 
ple have metallicities compiled from the literature by|Marbi^ 


et al. (20101. As noted by those authors, the metallicity data 


for many of these galaxies has been amassed over a number 
of years from a large number of disparate observations. As 
such it is challenging to put all the LVL metallicities into a 
common calibration framework. Indeed, Marble et al. ( 2010[ ) 
do not attempt to do so, simply noting that the non-uniform 
calibrations add an additional uncertainty across the sample 
of ~ 0.2 dex. We follow this approach, and do not attempt to 
standardise the LVL metallicity measurements. Instead, we 
apply an additional 0.2 dex of uncertainty to all LVL metal¬ 
licity measurements when fitting to the data. We stress that 
none of our conclusions are dependent on the inclusion of 
LVL - repeating our analysis without the LVL sample does 
not significantly change our results. 


2.3.2 Deriving H2 masses 

Molecular hydrogen masses can be calculated from the GO 
luminosity, Lqq, by assuming a GO/H 2 conversion factor 
aco- 


^(^2) — Q?coL/co(i-o) 


(3) 


There has been a vast amount of work over the last 
few decades dedicated to empirically measuring, and theo¬ 
retically modelling, aco ( Bolatto et al.|[2013 gives a recent 
review). While a number of physical factors can cause aco 
to vary, it is likely that the dominant factor is metallicity, 
with Qco increasing rapidly as metallicity decreases. In this 
work, we will present results by adopting a GO/H 2 conver¬ 
sion factor that varies with metallicity. Several metallicity- 
dependent CO/H 2 conversion factor prescriptions are given 
in a recent review by Bolatto et al. (20131. In particular. 


we will explore results derived using four recent metallicity- 
dependent conversion factors, de rived by [Wolfire et al.| 
(20101, Glover & Mac Low (2011 1, Narayanan et al. (20121, 
and |Feldmann et al.| ( |2012[ i. 

Uncertainties on molecular gas masses are dominated 
by systematic uncertainty concerning the GO/H 2 conversion 
factor. Uncertainties on GO fluxes, absolute flux calibration, 
and aperture corrections contribute, at most, ~ 20% uncer¬ 
tainty to the final value of M(H 2 ). The systematic uncer¬ 
tainty due to the CO/H 2 conversion factor is difflcult to 
estimate ( Saintonge et al.]|2011 estimate a combined uncer¬ 
tainty, including systematic uncertainty on the conversion 
factor, of 0.3 dex). We have attempted to mitigate system- 
atics due to the conversion factor by calculating our results 
for a range of metallicity-dependent factors. For the pur¬ 
poses of our analysis, however, we adopt a conservative total 
uncertainty of 0.3 dex on our molecular gas masses. 



log Stellar Mass 


Figure 3. The mass-metallicity relation for the galaxies in our 
sample. The grey shaded areas are taken from Tremonti et al. 
(2004), and show the area that contains 64% (light shaded area) 
and 90% (dark shaded area) of all SDSS galaxies. We have ad¬ 
justed the metallicity scaling for the Maiolino et al. (2008) calibra¬ 
tion used in this work. We have calculated molecular gas masses 
using the (metallicity-dependent) CO/H 2 conversion factor from 
Wolfire et al. (2010). The dashed line is the mass-metallicity fit 
to the sample of galaxies used in this work. Galaxies at high- 2 , 
which have greater uncertainty on their metallicity estimates, are 
shown with error bars. Open symbols show galaxies not detected 
in ^^CO, at the position of the 3a upper limits on their gas mass. 


3 RESULTS: MOLECULAR GAS AND THE 

MASS-METALLICITY RELATION 

We turn now to an examination of the effect of molecular gas 
on the scatter in the mass-metallicity relation . As discussed 
in the introduction, significant relationships have been dis¬ 
covered between the mass-metallicity relation and both SFR 
(the ‘FMR’; Mannucci et al.||2010 Lara-Lopez et aL||2010 1 
and Hi ( Bothwell et al.|2013 l. Any relationship with molec¬ 
ular gas has thus far proved elusive, due to a likely combina¬ 
tion of small available sample sizes, and potential degeneracy 
between metallicity and molecular gas masses (derived us¬ 
ing a conversion factor which is itself strongly dependent on 
metallicity). The combined sample assembled in this work 
represents a uniquely ideal dataset for examining any po¬ 
tential trends between metallicity and molecular gas mass, 
due to its dynamic range in both parameters. 

Figure shows the mass-metallicity relation for galax¬ 
ies in our combined sample. The original samples are dif¬ 
ferentiated with different symbols (shown in the inset key). 
The combined sample broadly follows the established 2 ~ 0 
mass-metallicity relation (the dark and light shaded regions 
in the background show the Icr and 2 a distributions in the 
mass-metallicity relation, as measured for SDSS galaxies by 
ITremonti et al.|2004| adjusted for the metallicity calibrations 
used in this work). As expected, higher -2 galaxies generally 
fall below the 2 ~ 0 mass-metallicity relation (we discuss this 
below). Most individual samples follow the mass-metallicity 
relation reasonably closely: several ALLSMOG galaxies, not 
detected in CO emission lie above the relation, while some 
LVL galaxies lie 0.2-0.5 dex below the relation. 

The symbol colour-coding in Fig.|^shows the molecular 
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Figure 4. The offset from the z = 0 mass-metallicity relation, plotting against molecular hydrogen mass {left panel), and total gas mass 
{right panel). In each plot, a linear fit is over-plotted. It can be seen that there is an inverse correlation between mass-metallicity offset 
and both H 2 mass and total gas mass (so galaxies with lower than average gas masses lie systematically above the mass-metallicity 
relation). 


gas mass, as derived using the Wolfire et al. (2010) CO/H 2 
conversion factor. A molecular gas ‘FMR’ effect would be 
apparent if galaxies at a constant stellar mass have an in¬ 
verse relation between molecular gas mass and metallicity 
Such a relation is difficult to discern by eye, due to both 
the uncertainty-driven scatter, and the the strong relation 
between molecular gas mass and stellar mass. To make any 
potential trends clearer, we examine the relation between 
molecular gas mass and the offset from the mass-metallicity 
relation (which we calculate by performing a quadratic fit to 
the mass-metallicity relation defined by our combined sam¬ 
ple). The quadratic fit used for calculating offsets is shown 
as a dashed line in Fig. (note that we have flattened the 
fit at upper end, to maximise at the ‘saturation metallicity’, 
12-|-log(0/H)=9.1). The fit is in good agreement with the 
fit found for a far larger SDSS sample by [Tremonti et al.| 
(20041, but is shifted slightly towards higher metallicities at 
stellar masses logM* > IO^'^Mq . 

Figure (left panel) shows the offset from the mass- 
metallicity relation, plotted against molecular gas mass. For 
legibility, we have only plotted error bars for the high -2 
galaxies, which have their metallicities measured using the 
[Nil]/Ha calibrator alone and as a result are more uncertain 
(due to the greater flux errors on the more distant objects, 
as well as the uncertainty due to the availability of just a 
single calibrator). 

We find that there is a shallow, but clear, inverse de¬ 
pendence between molecular gas mass and the offset from 
the mass-metallicity relation. Fitting linear functions to 
the relation, we find a slope of = —0.088 ± 0.021 (for the 


Wolfire et al. (20101 conversion factor; slopes derived using 


our other conversion factors are: Glover & Mac Low (20111 
= —0.078±0.019; Narayanan et al. ( 2012^ = — 0.098±0.020; 
and Feldmann et al. ) 20121 = —0.081±0.018). There is there¬ 
fore clear evidence for a molecular gas ‘FMR’ effect, whereby 


^ Note that we refer to the ‘FMR’ effect as being an inverse re¬ 
lation with metallicity at a given stellar mass; due to the low 
numbers of galaxies available with molecular gas measurements, 
we are unable to explore the potential redshift independence of 
our relations, another hallmark of the FMR. 


galaxies at a fixed stellar mass show an inverse dependence 
between their metallicity and their molecular gas mass. 

Figure (right panel) shows the offset from the mass- 
metallicity relation, plotted against total (i.e., H 2 -I-H 1 ) gas 
mass. For our high -2 galaxies, we have assumed that their 
ISM is entirely molecular (i.e., their total gas mass is equiv¬ 
alent to their molecular gas mass). As for the molecular gas 
above, we find that for each of four metallicity-dependent 
conversion factors there is a shallow inverse dependence be¬ 
tween total gas mass and the offset from the mass-metallicity 
relation. Given the dependence with molecular gas demon¬ 
strated above, and the dependence with atomic gas pre¬ 
sented by Bothwell et al. (2013), a dependence with total 
gas mass is of course unsurprising. Fitting linear functions 
the relation, we find a slope of = —0.076 ±0.022 (again, this 
is for the Wolfire et al. (20101 conversion factor; other slopes 
are:|GlovCT fc Mac Low ( |2011[ ) = —0.071 ±0.016; |Narayanan| 
et al. l 2012[ ) = —0.083±0.018; and Feldmann et al. (20121 = 


—0.077 ± 0.017). There is therefore also evidence for a total 
gas ‘FMR’ effect, whereby galaxies at a fixed stellar mass 
show an inverse dependence between their metallicity and 
their total (i.e., H 2 ±Hl) gas mass. 


3.1 Comparison to an analytical model 

In Fig. we have compared the right panel of Fig. to 
the prediction from the analytical ‘gas regulator’ model of 
Peng & Maiolino (2014). We have binned the model data 
(~ 5 X 10® model galaxies) in stellar mass bins of width 
0.5 dex, and colour-coded both the model tracks and the 
observational data by stellar mass. 

The model makes two main predictions. Firstly, the 
model predicts that within a stellar mass bin, the offset from 
the mass-metallicity relation does increase (in the negative 
direction, i.e., below the mass-metallicity relation) with in¬ 
creasing total gas mass. This is broadly the same trend as we 
have found for our observational data above. However, the 
model also demonstrates that the dependency of the mass- 
metallicity offset is far steeper within any individual stellar 
mass bin than it would be for the sample as a whole (due 
to the strong correlation between stellar mass and total gas 
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Figure 5. The offset from the z = 0 mass-metallicity relation, plotted against total gas mass, for (left panel) our observed sample, and 
(right panel) calculated using the Peng &: Maiolino (2014) gas regulator model. Model data are shown in bins of stellar mass of 0.5 dex 
width. The model predicts a strong inverse dependence between mass-metallicity offset within each stellar mass bin. Also shown is the 
predicted ‘quenching sequence’ (dashed line, with associated standard deviation indicated with dotted lines). 
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mass). However, such a separation in stellar mass is not im¬ 
mediately clear in the observational data - there seems to 
be sufficient scatter in the gas mass-stellar mass correlation 
that any potentially clear trends are obscured. The model 
does demonstrate why, when considering the the sample as a 
whole, we only see a shallow dependence between metallicity 
offset and gas mass. 


4 RESULTS: PRINCIPLE COMPONENT 

ANALYSIS 

The analysis contained within the previous sections suffers 
from two potential weaknesses. Firstly, as previously dis¬ 
cussed, the correlation between gas mass and stellar mass 
serves to erode any clear trends between mass-metallicity 
offset and gas mass; other than showing that the analyti¬ 
cal model is capable of reproducing the observed data, it 
is challenging to quantify any trends. Secondly, the tech¬ 
nique suffered from being parametric - that is, dependent on 
the assumption that some distribution exists between stellar 
mass and metallicity, and paramerising data points in terms 
of their offset from a quadratic fit to that relation. In our 
case, the specific results obtained depend strongly on the 
form of the fit to the mass-metallicity relation. Ideally, we 
would like to frame the analysis in terms of a non-parametric 
fitting technique, which simply examines the correlations in 
the data without presupposing any specific relationship be¬ 
tween them. One snch technique is Principle Component 
Analysis (PCA). 

Principle Component Analysis (PCA) is a linear trans¬ 
formation technique, by which a set of physical variables 
are converted into a set of orthogonal, linearly uncorrelated 
components (called ‘Principle Components’) defined by a 
set of eigenvectors. The conversion is chosen so to ensure 
that the maximum possible variance in the data is contained 
within the first component, and each subsequent component 
contains as much variance as possible (with the constraint 
that each component is orientated orthogonally to every 
preceding component). In practice, PCA performs a coor¬ 
dinate transformation revealing the optimum ‘projection’ of 


a dataset, and revealing which parameters are responsible 
for the variance in the sample. PCA is particularly useful 
for revealing any possible ‘dimensionality reduction’ - for 
example, revealing that a distribution forms a 2D ‘plane’ in 
3D parameter space (as [Hunt et al.|[2012| found for stellar 


mass, SFR and metallicity). PCA was used by both Lara- 


Lopez et al. (20101 and Hunt et al. (20121 to examine the 


secondary dependence of SFR on the mass-metallicity rela¬ 
tion. 


Here, we perform PCA on our data, examining the stel¬ 
lar mass, metallicity, and a third parameter of interest. In 
turn, we will examine the molecular gas mass, the total gas 
mass, the star formation rate, and the star formation effi¬ 
ciency. We first normalise each parameter to the mean value 
for our combined sample: 


log (M.)'’'^^ = log (M.) - 10.13 
12 -f log(0/H)^‘^^ = 12 + log(0/H) - 8.95 
log (SFE)^'^^ = log (SFE) + 8.68 
log (SFR)'’^^ = log(SFR) - 0.39 
log (Mh 2 )'’°^ = log (Mh 2 ) - 9.07 

log (Mg,,)P‘^^ = log (Mgas) - 9.71 

We have accounted for uncertainty by adopting a 
Monte-Carlo approach, performing 10® PCA iterations for 
each combination of parameters. For each iteration, each 
galaxy has the value of its physical parameters randomly 
perturbed by an amount following the respective error on 
each parameter. At the end of the 10® iterations, we take 
the mean and standard deviation of the resulting Eigen¬ 
vector distribution to be our Hnal Eigenvector values (and 
uncertainties thereon). In order to ensure we are not unduly 
influenced by outliers, we also perform sample bootstrap¬ 
ping during the Monte-Carlo: for each iteration, we ran¬ 
domly sample (with replacement) our complete sample of 
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PC1 = 0.578 (M*) - 0.807 (MH2) + 0.118 (Z) 
PC2 = - 0.760 (M*) - 0.585 (MH2) - 0.280 (Z) 
PCS = 0.295 (M*) - 0.072 (MH2) - 0.952 (Z) 


Variation in PCI = ( 79.7 +/- 0.6) % 
Variation in PC1+PC2 = ( 97.0 +/- 0.3) % 
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Figure 6. Three projections of the 3-space defined by the principle components derived in 0 The inset legends above the plots give 
the definitions of the principle components, in terms of stellar mass (M*), molecular gas mass (MH2), and metallicity (Z). 


galaxies, generating in each case a sample size equal to our 
original dataset. 

It must be noted that a potential weakness in the ap¬ 
plication of PCA to our data is that PCA can only describe 
datasets in terms of linear relationships between parameters 
- more complex distributions of data cannot be described 
in terms of a simple set of orthogonal eigenvectors. While 
linear relationships (such as the SFR-M* ‘main sequence’) 
are ideal for modelling using PCA, applying the technique 
to non-linear correlations (such as the mass-metallicity re¬ 
lation) will, be definition, be somewhat inaccurate. In prac¬ 
tice, this will have the effect of increasing the apparent scat¬ 
ter around the component vectors. Given the relatively low 
numbers of galaxies in our sample, it is likely that the un¬ 
certainty added by forcing linear relations onto the data is 
not larger than the uncertainty inherent in the distribution 
of the data (as revealed by our Monte-Carlo approach). We 
therefore continue to perform PCA on our data, but cau¬ 
tion the reader that some uncertainty may be added due to 
non-linearity in some underlying correlations. 

The results are discussed below. Fig.j^shows the results 
of performing PCA on the parameters stellar mass, metal¬ 
licity, and molecular gas mass. Appendix Figures and[^ 
show the results when using total gas mass, star formation 
rate, and star formation efficiency (respectively) as the third 


parameter. Tablej^lists the values of the Eigenvectors of the 
components in each instance. 

In each of the four cases, the PCA results do indeed 
suggest that the galaxy form a 2D plane in the parameter 
space, with just ~ 1.5% —3% of the total variation occurring 
in the third principle component. As noted by [Hunt et al.| 
(20121, with very little of the sample variance contained 
within the third principle component, it is a useful tool 
for examining the optimal relationship between the three 
parameters of interest. In each of our cases, the third pa¬ 
rameter is highly dominated by metallicity, with only minor 
contributions from the other parameters. We can therefore 
‘solve’ for metallicity, by setting the third principle compo¬ 
nent equal to zero (which is valid at the ~ 1.5% — 3% level), 
allowing us to examine the effect of each parameter (stellar 
mass, molecular gas mass, total gas mass, SFR, and SFE) 
on metallicity. 


Molecular gas mass: 

The PCA results for molecular gas are shown in Fig. 

~ 80% of the variance is contained within the first compo¬ 
nent, and ~ 97% of the variance is contained within the first 
two. The first component is defined mainly in terms of stel¬ 
lar mass and molecular gas mass (molecular gas mass being 
the dominant contributor), with a small contribution from 
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PC1 = 0.441 (M*) - 0.893 (Mgas) + 0.085 (Z) 
PC2 = - 0.858 (M*) - 0.448 (Mgas) - 0.249 (Z) 
PCS = 0.261 (M*) - 0.037 (Mgas) - 0.964 (Z) 


Variation in PCI = ( 72.9 +/- 0.6) % 
Variation in PC1+PC2 = ( 97.3 +/- 0.2) % 
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Figure 7. Three projections of the 3-space defined by the principle components derived in 0 The inset legends above the plots give 
the definitions of the principle components, in terms of stellar mass (M*), total gas mass (Mgas), and metallicity (Z). 


metallicity. The second component is similarly defined, with 
the contribution from metallicity being a little greater. The 
third component, which accounts for just 2.8% of the sample 
variation, is dominated by metallicity (telling us that that 
metallicity has the minority contribution to the total sample 
variance). 

We now set the third principle component to zero. Here, 
the third principle component is defined as 

PCS = 0.295 (M*)'’'^^ - 0.072 (Mhz)^^^ 

- 0.952 12 + log(0/H)^'^^ (4) 

Setting this to zero implies that: 

12 + log(0/H) = 0.31(log M.) - 0.08(log Mhz) + 6.53 (5) 

That is, the metallicity is primarily determined by the 
stellar mass (i.e., the mass-metallicity relation), with a sec¬ 
ondary dependence on the molecular gas mass. This sec¬ 
ondary dependence on molecular gas mass is approximately 
22% times as strong as the primary dependence on stellar 
mass. 

Total gas mass: 

The PCA results for total gas mass are shown in Fig. 


~ 73% of the variance is contained within the first compo¬ 
nent, and ~ 98% of the variance is contained within the first 
two. The first component is defined mainly in terms of stel¬ 
lar mass and total gas mass, with a small contribution from 
metallicity. And again, the second component is similarly 
defined, with the contribution from metallicity being a little 
greater. The third component, which here accounts for just 
2.5% of the sample variation, is dominated by metallicity. 

Again, we set the third principle component to zero. 
Here, the third principle component is defined as 

PC3 = 0.261 (M*)'’^^ - 0.038 (Mga.s)^'^^ 

- 0.964 12 -b log(0/H)^°^ (6) 

Setting this to zero implies that: 

12 -b log(0/H) = 0.27(log M*) - 0.04(log Mg,,) -b 6.60 (7) 

That is, the metallicity is primarily determined by the 
stellar mass, with a secondary dependence on the total gas 
mass. This secondary dependence on total gas mass is ap¬ 
proximately 13% times as strong as the primary dependence 
on stellar mass. This shows that the metallicity dependence 
with the molecular gas content (discussed in the previous 
sub-section) is stronger than with the total gas content. 
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PC1 = 0.513 (M*) - 0.852 (SFR) + 0.100 (Z) 
PC2 = - 0.816 (M*) - 0.521 (SFR) - 0.249 (Z) 
PCS = 0.265 (M*) - 0.045 (SFR) - 0.962 (Z) 


Variation in PCI = ( 72.6 +/- 0.4) % 
Variation in PC1+PC2 = ( 97.0 +/- 0.3) % 
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Figure 8. Three projections of the 3-space defined by the principle components derived in 0 The inset legends above the plots give 
the definitions of the principle components, in terms of stellar mass (M*), star formation rate (SFR), and metallicity (Z). 


Star formation rate: 


The PCA results for SFR are shown in Fig. ~ 73% 
of the variance is contained within the first component, and 
~ 97% of the variance is contained within the first two. 
The first component is defined mainly in terms of stellar 
mass and SFR, with a small contribution from metallicity. 
And again, the second component is similarly defined, with 
the contribution from metallicity being a little greater. The 
third component, which here accounts for just 2.8% of the 
sample variation, is dominated by metallicity. These results 


are consistent with the work done by Hunt et al. (20121, 


who used PCA to examine correlations between stellar mass, 
SFR and metallicity for a sample of metal-poor starbursts. 

Setting the third principle component to zero. Here, the 
third principle component is defined as 


PCS = 0.265 (M*)'’^^ - 0.045 (SFR)'’^^ 

- 0.963 12 + log(0/H)^°^ (8) 

Setting this to zero implies that: 


12 + log(0/H) = 0.28 (log M*) - 0.05(log SFR) -f 6.14 (9) 


stellar mass, with a secondary dependence on the SFR. This 
secondary dependence on SFR is approximately 16% times 
as strong as the primary dependence on stellar mass. As 
before, this dependence on the SFR is weaker than the de¬ 
pendence on the molecular gas content. 

Star formation efficiency: 

The PCA results for SFE are shown in Fig. ~ 66% 
of the variance is contained within the first component, and 
~ 97% of the variance is contained within the first two. The 
first component is defined mainly in terms of stellar mass 
and SFE, with a slightly smaller contribution from metal¬ 
licity. The second component is similarly defined, with the 
contribution from metallicity being a somewhat less. The 
third component, which here accounts for 3.0% of the sam¬ 
ple variation, is dominated by metallicity, with only a very 
minor contribution from SFE. 

Again, we set the third principle component to zero. 
Here, the third principle component is defined as 


PC3 = 0.233 (M*)'”'^^ - 0.007 (SFE)^'^^ 

- 0.972 12 -b log(0/H)'’'^^ (10) 


That is, the metallicity is primarily determined by the 


Setting this to zero implies that: 
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PC1 = 0.904 (M*) - 0.368 (SFE) + 0.214 (Z) 
PC2 = - 0.357 (M*) - 0.929 (SFE) - 0.092 (Z) 
PC3 = 0.233 (M*) - 0.007 (SFE) - 0.972 (Z) 


Variation in PCI = ( 64.7 +/- 0.7) % 
Variation in PC1+PC2 = ( 95.1 +/- 0.4) % 
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Figure 9. Three projections of the 3-space defined by the principle components derived in ij^ The inset legends above the plots give the 
definitions of the principle components, in terms of stellar mass (M*), star formation efficiency (SFE), and metallicity (Z). Metallicity 
has a weaker effect on the distribution, which can be mainly defined in terms of M* and SFE. 


12 + log(0/H) = 0.24(log M.) - 0.008(log SFE) + 6.45 

( 11 ) 

That is, the metallicity is primarily determined by the 
stellar mass, with negligible secondary dependence on the 
SFE: the secondary dependence on SFE is approximately 3% 
times as strong as the primary dependence on stellar mass. 
This clearly indicates that there is no correlation between 
SFE and metallicity, once the dependence on stellar mass 
is taken into account. This will be discussed further in the 
next section. 


5 DISCUSSION 

5.1 The effect of a changing CO/H 2 conversion 
factor 


Thus far, all principle component analysis results have been 
calculated using the metallicity-dependent CO/H 2 conver¬ 
sion factor presented by Wolfire et al. (20101. It is, however, 
important to investigate the extent to which these results 
are robust to a range of conversion factor prescriptions. 

We perform the same PCA steps for three further metal¬ 


licity dependent CO/H 2 conversion factors: Glover & Mac 


Low (20111, Feldmann et al. (20121, and Narayanan et al. 


(20121. In the interest of saving space, we do not repeat 
the above discussion for each case. Nevertheless, each con¬ 
version factor produces highly similar behaviour: examining 
each secondary-dependence in turn (molecular gas mass, to¬ 
tal gas mass, SFR, and SFE), we find that the distribution 
of data does form a 2D plane in 3D parameter space, with 
between 2% — 5% of the sample variation being contained 
within the first two principle components. 

We focus on the differences produced by the final step in 
the analysis - that by setting the third principle component 
to zero, the metallicity can be defined in terms of the stellar 
mass, with a secondary modifying term coming from one of 
the parameters listed above (in the original formulation of 
the FMR by [Mannucci et al.]|2010| and |Lara-L6pez et J)] 


2010 SFR was used for this third parameter). 


As above, we parameterise the strength of the effect of 
the third parameter (either M(H 2 ), M(gas), SFR, or SFE) on 
the metallicity by setting the third principle component to 
zero, and finding the resultant expression for the metallicity: 


12 + log(0/H) =log(M*)-/rlog(X), (12) 

where X is, in turn, M(H 2 ), M(gas), SFR, or SFE. Fig. 
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Component 

M* 

M(H2) 

12-tlog(0/H) 

PC 1 

PC 2 

PC 3 

±0.578 ± 0.040 
-0.760 ± 0.033 
±0.295 ± 0.027 

-0.807 ± 0.032 
-0.585 ±0.043 
-0.072 ±0.021 

±0.118 ±0.023 
-0.280 ±0.031 
-0.952 ±0.008 

Component 

M* 

M(gas) 

12±log(0/H) 

PC 1 

PC 2 

PC 3 

-1-0.441 ± 0.081 
-0.858 ±0.045 
±0.261 ± 0.024 

-0.893 ±0.041 
-0.448 ± 0.083 
-0.037 ± 0.015 

±0.085 ± 0.023 
-0.249 ±0.022 
-0.964 ± 0.007 

Component 

M* 

SFR 

12±log(0/H) 

PC 1 

PC 3 

PC 3 

±0.513 ±0.056 
-0.816 ±0.036 
±0.265 ± 0.021 

-0.852 ±0.036 
-0.521 ± 0.059 
-0.045 ±0.014 

±0.100 ±0.027 
-0.249 ± 0.021 

-0.962 ± 0.005 

Component 

M* 

SEE 

12±log(0/H) 

PC 1 

PC 2 

PC 3 

±0.904 ± 0.042 
-0.357 ±0.113 
±0.233 ±0.025 

-0.368 ±0.114 
-0.929 ± 0.044 
-0.007 ± 0.004 

±0.214 ± 0.023 
-0.092 ±0.032 
-0.972 ±0.006 


Table 1. Eigenvectors of the three principle components, for each 
of our parameter sets (stellar mass, metallicity, and a third pa¬ 
rameter of interest - respectively, molecular gas mass, total gas 
mass, SFR, and SEE. Errors have been obtained via a Monte 
Carlo bootstrapping method, as described in the text. 


|10| shows the values of we find for each parameter. All 
values of are < 1, implying that in all cases the dominant 
driver of metallicity variations is the stellar mass. For each 
of the four CO/H 2 conversion factors, we find that gas con¬ 
tent is the most influential ‘secondary parameter’. For three 


out of four of the conversion factor prescriptions (Wolflre 


et al.||20T0| [Glover fc Mac Low||20lH and [Narayanan et al 


20121, it is the molecular gas that has the most influence on 


metallicity (after stellar mass). For the remaining conver¬ 
sion factor prescription - Feldmann et al. (20121 - it is the 
total gas content that is the strongest secondary parameter 
(though it is comparable to the effect of molecular gas within 
the uncertainties on jj,, which are typically ±0.05). But, in¬ 
dependent of the choice of conversion factor, we find that 
the strongest ‘FMR’ effect is to be found with gas content. 


5.2 Comparison to the SFR-FMR 


The strength of SFR as a secondary parameter is somewhat 
less than gas content, with (j-sfr ~ 0.17 (small variations 
in fj-sFR shown in Fig. are due to randomness introduced 
by the bootstrapping technique). We note that this value is 
significantly smaller than the value of ^ found by|Mannucci[ 


et al. (20101: they express metallicity as a function of stellar 


mass and SFR: 


12 ± log(0/H) = log (M*) - 0.32 log(SFR) (13) 

We therefore find the size of the effect of SFR on metal- 
approximately 50% of that found by |Mannucci[ 
). This discrepancy is far larger than can be ac¬ 
counted for by sample uncertainty; our bootstrapping anal- 


hcity to be 


et al. (2010 


ysis estimates a combined uncertainty on the influence of 
SFR of ~ 30% (AtsFR = 0.17 ±0.05). 

One potential resolution to the conflict lies in the 
method used to calculate the SFRs. Mannucci et al. (20101 
took their sample from SDSS, using the SFRs and strong line 
fluxes as provided by the SDSS spectroscopic fibre. SDSS 
calculates its pipeline SFRs by fitting to optical strong lines 
(with the Hq having the most influence). As such (a) previ¬ 
ously, SFRs and metallicities were calculated from matched- 
aperture observations, and (b) the same optical strong lines 
were used to calculate both SFRs and metallicities, which 
could introduce a co-dependency. Both of these effects could 
result in the correlation between SFR and metallicity being 
artificially tightened. It has been shown that when calculat¬ 
ing SFRs by adopting integrated SFRs, corrected for fibre- 
aperture effects, the SFR-FMR effect weakens considerably 
(Bothwell et al.|[2013 Zahid et al.|2014|. 


In this work, we have used both galaxies with SDSS- 
SFRs (galaxies taken from the ALLSMOG and COLD GASS 
surveys), and galaxies with SFRs calculated with other 
methods (we used a combination of IR and UV fluxes for 
the HRS and LVL samples, the z ~ 2 disks have SFRs de¬ 
rived from fitting to a combination of UV, optical contin¬ 
uum, Ha, and 24^m continuum fluxes, and the z ~ 2 SMGs 
have SFRs derived using IR fluxes, via the far IR-radio cor¬ 
relation) . These disparate methods, while all being effective 
probes of the total star formation rate, will clearly not intro¬ 
duce the same co-dependency as taking all values from the 
same SDSS spectra. This weaken any correlations between 
SFR and metallicity which were previously spuriously en¬ 
hanced. 


5.3 Comparison with an analytical model 

A key prediction of the model we discuss in §3.2.1 is the exis¬ 
tence of a ‘quenching sequence’, populated by galaxies which 
are elevated above the mass-metallicity relation. Within the 
framework of the model, the only mechanism able to ele¬ 
vate galaxies significantly above the mass-metallicity rela¬ 
tion is ‘strangulation’ - the cessation of gas infall - which 
leads firstly to the build up of metals due to star forma¬ 
tion (with a lack of metal-poor inflowing gas acting to di¬ 
lute the ISM), and ultimately to the shutting down of star 
formation (‘quenching’). Recently, Peng et al. (2015) have 
shown that the comparison of stellar metallicities between 
passive and star forming galaxies strongly support strangu¬ 
lation as the main quenching mechanism for galaxies with 
mass < IO^^Mq. Being an emission line-selected sample, 
none of our sample galaxies are currently fully quenched. 
The model predicts, however, that the ~ 7 galaxies elevated 
above our relation are essentially analogues of ‘green val¬ 
ley’ galaxies, identified via their metallicity and gas mass, 
which have recently undergone strangulation and are in the 
processes of becoming quenched and dead. 

It is important to investigate whether the model is doing 
a good job of explaining the scatter in the data. Randomness 
in the model is introduced by allowing the gas inflow rate 
to vary by a factor of two - this introduces scatter in both 
metallicity and gas content. If it can be seen that the scatter 
in the observed data is greater than that contained in the 
model (allowing for observational uncertainties), it is likely 
that real gas inflow rates vary by a factor of > 2. We compare 
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Figure 10. The strength of the influence on the mass-metallicity 
relation of molecular gas mass (green line), total gas mass (blue 
line), SFR (red line), and SFE (black line). The strength of the 
influence is quantified as /r — see Eq. 12. Results for four recent 
metallicity-dependent CO/H 2 conversion factors are shown. In 
each case, the strongest effect is found to be with gas content 
(either molecular or total). The smallest influence on the mass- 
metallicity relation is found with star formation efficiency. 


the model to the data as follows. For each observed galaxy, 
we identify all model galaxies with the same stellar mass 
and metallicity (within bins of 0.05 dex). We then compare 
the total gas mass of the real galaxy to the mean gas mass of 
the (mass & metallicity matched) model galaxies. We find a 
mean discrepancy between the model and observed gas mass 
of 0.49 dex. This must be compared, though, to the variance 
in gas mass contained within the model itself. Within each 
0.05 dex bin of stellar mass and metallicity, the standard 
deviation of model gas masses is 0.41 dex. That is, the mean 
discrepancy between the analytical model and the data is 
only slightly more than the variance contained within the 
model itself. Allowing for typical observational uncertainties 
on the gas mass the model provides a good match to the 
data, suggesting that the assumptions of the model (i.e., that 
gas inflow rate varies by up to a factor of 2) are sufficient to 
explain the observed scatter. 


tion efficiency for our sample would result in a good, yet 
spurious, correlation. 


At first glance, a lack of correlation between metallicity 
and star formation efficiency appears somewhat surprising: 
metal lines act as an efficient cooling mechanism for the ISM, 
and given a lack of metals the ISM cannot effectively cool, 
fragment into molecular clouds, and form stars. This was the 
conclusion drawn by Shi et al. (20141, who found extremely 
low star formation efficiencies in two local metal-poor dwarf 
galaxies (these low SFEs result from a large reservoir of CO- 
dark molecular gas, which |Shi et aT]|2014| trace using dust 
emission). However, the low metallicity regime in which the 


Shi et al. (20141 study was carried out may differ signih- 


cantly from the ISM conditions within the galaxies analysed 
in this work, and it is unclear if the results are applica¬ 
ble; many models predicting reduced SFE due to inefficient 


cooling in low metallicity environments (i.e., Krumholz et al. 


20091 deal with galaxies at such low metallicities that the 


ISM is entirely atomic. 


There are few explorations of the connection between 
SFE and metallicity in the molecular regime. [Dib' et al.| 
(20111 present a semi-analytic model of protocluster clumps 
which predicts an increase in star formation efficiency as 
metallicity decreases, primarily due to the weaker stellar 
winds at low metallicity (which results in a less thorough 
evacuation of gas in low-metallicity star-forming clumps). 
Observationally, Dib et al. (20111 find a trend supporting 
this, which is, however, based on just three galaxies (the 
SMC, M33, and NGC6822). 


To date, therefore, there has been no clear observational 
consensus as to the effect of metallicity on star formation ef- 
hciency, with previous studies having very small sample sizes 
(2-3 galaxies) and reaching conflicting results. The results 
presented in this work, which find a negligible connection 
between star formation efficiency and metallicity, represent 
the first investigation of this issue with a significant sample 
size. 


5.4 The connection between metallicity and star 
formation efficiency 

In all cases, by far the weakest ‘FMR’ effect is found with 
the star formation efficiency. Across all CO/H 2 conversion 
factors, the SFE and metallicity are only very weakly con¬ 
nected. In other words, the Schmidt-Kennicutt relation does 
not depend on metallicity. To date, the relationship between 
metallicity and star formation efficiency has been challeng¬ 
ing to examine, due to the strong correlation between metal¬ 
licity and stellar mass. The molecular gas consumption time 


{th 2 = 1/SEE) is observed to vary with stellar mass (Sain- 
tonge et al.|201l[ Bothwell et al. 20141, a phenomenon which 


Saintonge et al. (20111 attribute to the presence of increas¬ 


ingly ‘bursty’ star formation histories towards lower halo 
masses. It is due to the strong co-variance between metal¬ 
licity and stellar mass that a technique like Principle Com¬ 
ponent Analysis is required, in order to break the mass- 
metallicity degeneracy and isolate the effects of individual 
parameters. Simply plotting metallicity against star forma- 


There are also many physical effects potentially govern¬ 
ing the effect of metallicity on star formation efficiency (gas 
cooling and stellar wind efficiency being just two), and no 
clear theoretical consensus on which effects will dominate. 
Phenomenologically, the weakness of the connection between 
SFE and metallicity follows simply from the fact that metal¬ 
licity varies with molecular gas content and SFR in similar 
ways; a secondary dependence in the mass-metallicity rela¬ 
tion is observed with both molecular gas and SFR (though 
molecular gas has the stronger effect). As this dependence is 
in the same direction for each parameter - at a given stellar 
mass, galaxies with either increased SFR or increased M(H 2 ) 
have lower metallicity - when considering the ratio of these 
two quantities (the SFE), the effect nearly cancels out. As a 
result, SFE and metallicity are only very weakly connected 
in the samples considered in this work. Physically, it may 
be that a combination of metallicity-dependent processes 
(e.g. weak stellar winds, inefficient gas cooling) conspire to 
leave no connection between SEE and metallicity. Such a 
conclusion, however, would require theoretical investigation 
beyond the scope of this work. 
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6 CONCLUSIONS 

In this work, we have presented an analysis of the connec¬ 
tion between the gas content (both molecular and total) and 
metallicity of galaxies in samples from 0 < z < 2. We have 
explored the effect of various metallicity-dependent CO/H 2 
conversion factors on the Schmidt-Kennicutt relation. Con¬ 
versely, we have also investigated the effect of gas content 
on the mass-metallicity relation, the most commonly used 
metallicity scaling relation. Finally, we use Principle Com¬ 
ponent Analysis to statistically analyse these parameters 
within a non-parametric framework. Our main conclusions 
are as follows: 

• We demonstrate the existence of a systematic ‘Fun¬ 
damental Metallicity Relation’ effect between molecular 
gas mass, stellar mass, and metallicity, by which at a fixed 
stellar mass, galaxies with higher H 2 mass lie systematically 
below the mass-metallicity relation. 

• We use Principle Component Analysis to explore the 
correlations between stellar mass, metallicity, and a third 
parameter (respectively, molecular gas content, total gas 
content, SFR, and SFE) in a non-parametric way. We find 
strong evidence for secondary dependences with both gas 
content and SFR, with the effect of gas content being 
stronger independent of the choice of CO/H 2 conversion 
factor. We conclude that gas content is the driver of 
the ‘Fundamental Metallicity Relation’, with the (more 
commonly used) SFR-FMR being simply a consequence 
due to the Schmidt-Kennicutt relation. 

• We find a very negligible connection between gas- 
phase metallicity and star formation efficiency; that is, the 
Schmidt-Kennicutt relation has no dependence on metallic¬ 
ity. This result stands in opposition to some recent findings, 
which were conducted with far smaller sample sizes than 
used in this work. 
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log M* 
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12+log(0/H) 
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8.89 

1.70 
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8.34 
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0.314 

<8.08 

ALLSMOG 

8.96 

8.65 

1.50 

<8.04 

ALLSMOG 

8.54 

8.55 

0.0357 

<8.50 

ALLSMOG 

9.06 

8.79 

12.0 

<7.95 

ALLSMOG 

10.17 

8.95 

2.11 

8.87 

COLD GASS 

10.05 

8.93 

3.47 

8.87 

COLD GASS 

10.08 

8.99 

2.19 

8.91 

COLD GASS 

10.81 

9.15 

1.35 

9.14 

COLD GASS 

10.74 

9.17 

4.41 

9.13 

COLD GASS 

10.67 

9.07 

7.40 

9.73 

COLD GASS 

10.74 

9.07 

9.33 

9.85 

COLD GASS 

10.90 

9.10 

6.59 

9.81 

COLD GASS 

10.26 

9.08 

0.910 

8.36 

COLD GASS 

10.41 

9.15 

2.48 

8.80 

COLD GASS 

10.55 

9.13 

2.41 

9.14 

COLD GASS 

10.37 

9.03 

3.40 

9.21 

COLD GASS 

10.48 

9.10 

1.13 

9.03 

COLD GASS 

10.10 

9.09 

1.38 

8.96 

COLD GASS 

10.53 

9.13 

26.8 

9.30 

COLD GASS 

10.07 

8.88 

2.18 

8.85 

COLD GASS 


Table 2. Physical parameters (stellar mass, metallicity, SFR, and M(H 2 ) derived using the Wolfire et ah 
factor) for the galaxies in this work. 


(2010) CO/H 2 conversion 


© 2014 RAS, MNRAS 000, ??-?? 



18 M. S. Bothwell et al. 


log M* 
[M©] 

12+log(0/H) 

SFR 

[Mo/yr] 

M(H2) 

[M©] 

Source 

10.11 

9.08 

2.42 

9.15 

COLD GASS 

10.55 

9.06 

9.75 

9.59 

COLD GASS 

10.37 

9.04 

8.57 

9.46 

COLD GASS 

10.06 

9.00 

1.94 

8.97 

COLD GASS 

10.11 

9.08 

2.53 

8.93 

COLD GASS 

10.20 

9.00 

0.810 

8.90 

COLD GASS 

10.38 

9.10 

2.89 

9.35 

COLD GASS 

10.01 

9.07 

1.28 

8.68 

COLD GASS 

10.32 

8.87 

10.0 

9.58 

COLD GASS 

10.98 

9.07 

3.75 

9.84 

COLD GASS 

10.27 

9.02 

1.00 

9.17 

COLD GASS 

10.64 

9.09 

3.86 

9.52 

COLD GASS 

10.03 

8.95 

1.24 

9.49 

COLD GASS 

10.59 

9.05 

7.31 

9.63 

COLD GASS 

10.24 

9.10 

0.110 

8.75 

COLD GASS 

10.02 

8.98 

0.990 

8.60 

COLD GASS 

10.39 

9.04 

3.38 

9.29 

COLD GASS 

10.44 

9.00 

16.0 

9.51 

COLD GASS 

10.80 

9.17 

2.72 

9.18 

COLD GASS 

10.29 

9.13 

1.45 

8.80 

COLD GASS 

10.98 

9.15 

0.620 

8.91 

COLD GASS 

10.70 

9.14 

3.61 

9.53 

COLD GASS 

10.09 

9.05 

5.72 

9.23 

COLD GASS 

10.17 

9.07 

2.41 

8.81 

COLD GASS 

10.77 

9.07 

5.51 

9.68 

COLD GASS 

10.91 

9.05 

7.11 

9.39 

COLD GASS 

10.60 

9.18 

3.46 

9.07 

COLD GASS 

10.49 

9.05 

13.7 

9.72 

COLD GASS 

10.15 

9.03 

2.11 

8.86 

COLD GASS 

10.12 

9.06 

2.03 

9.02 

COLD GASS 

10.77 

9.04 

6.84 

9.76 

COLD GASS 

10.28 

9.00 

9.40 

9.63 

COLD GASS 

11.03 

9.08 

4.80 

9.61 

COLD GASS 

10.85 

8.97 

3.50 

9.61 

COLD GASS 

10.95 

9.19 

3.47 

9.63 

COLD GASS 

10.15 

9.05 

2.60 

9.17 

COLD GASS 

10.99 

8.92 

0.410 

9.97 

COLD GASS 

10.31 

9.04 

4.00 

9.28 

COLD GASS 

10.18 

8.96 

9.73 

9.42 

COLD GASS 

10.44 

9.15 

9.92 

9.49 

COLD GASS 

11.33 

9.14 

1.92 

9.35 

COLD GASS 

10.91 

9.08 

2.77 

9.74 

COLD GASS 

10.41 

9.15 

3.22 

9.10 

COLD GASS 

10.18 

8.92 

1.01 

9.02 

COLD GASS 

10.76 

9.07 

5.74 

9.56 

COLD GASS 

10.42 

9.20 

4.50 

9.34 

COLD GASS 

10.77 

9.05 

4.72 

9.73 

COLD GASS 

10.87 

9.08 

5.43 

9.53 

COLD GASS 

10.57 

9.10 

1.93 

8.96 

COLD GASS 

10.54 

9.16 

1.46 

8.87 

COLD GASS 

11.03 

9.18 

1.47 

8.91 

COLD GASS 

10.46 

9.09 

0.240 

8.51 

COLD GASS 

10.28 

9.13 

2.41 

8.77 

COLD GASS 

10.04 

9.06 

2.40 

9.04 

COLD GASS 

10.02 

8.99 

0.350 

8.43 

COLD GASS 

10.05 

9.03 

1.24 

8.90 

COLD GASS 

10.56 

9.16 

1.72 

8.67 

COLD GASS 

10.64 

9.14 

7.00 

9.69 

COLD GASS 

10.10 

8.98 

2.14 

9.27 

COLD GASS 

10.07 

9.10 

1.43 

8.59 

COLD GASS 

10.13 

9.06 

1.57 

8.98 

COLD GASS 

10.05 

9.03 

1.68 

8.90 

COLD GASS 

10.60 

9.02 

3.44 

9.51 

COLD GASS 

10.25 

9.14 

1.45 

8.67 

COLD GASS 


Table 2. (cont) Physical parameters (stellar mass, metallicity, SFR, and M(H 2 ) derived using the Wolfire et al. (2010) CO/H 2 conversion 
factor) for the galaxies in this work. © 2014 RAS, MNRAS 000, ??-?? 
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log M, 
[Mo] 

12+log(0/H) 

SFR 

[Mo/yr] 

M(H2) 

[Mq] 

Source 

10.09 

8.96 

0.900 

<8.88 

COLD GASS 

11.18 

9.14 

2.25 

<9.17 

COLD GASS 

10.24 

9.04 

0.350 

<8.93 

COLD GASS 

10.01 

9.20 

2.17 

<9.10 

COLD GASS 

10.03 

9.12 

0.0600 

<8.83 

COLD GASS 

10.09 

8.94 

0.000 

<8.62 

COLD GASS 

11.01 

9.12 

0.520 

<9.16 

COLD GASS 

10.89 

9.17 

0.0400 

<9.09 

COLD GASS 

11.13 

9.19 

0.430 

<9.25 

COLD GASS 

10.79 

9.14 

1.19 

<9.12 

COLD GASS 

11.08 

9.15 

0.790 

<9.18 

COLD GASS 

10.11 

8.95 

0.0600 

<8.65 

COLD GASS 

9.46 

8.81 

3.04 

8.89 

Herschel Reference Survey 

9.69 

8.94 

6.00 

8.67 

Herschel Reference Survey 

10.2 

9.02 

13.8 

9.34 

Herschel Reference Survey 

10.2 

8.96 

16.5 

9.23 

Herschel Reference Survey 

9.77 

8.76 

8.34 

9.16 

Herschel Reference Survey 

10.0 

9.01 

5.02 

9.10 

Herschel Reference Survey 

9.42 

8.83 

3.44 

8.45 

Herschel Reference Survey 

9.25 

8.68 

1.05 

8.62 

Herschel Reference Survey 

9.28 

8.85 

1.80 

8.47 

Herschel Reference Survey 

9.20 

8.79 

1.10 

8.34 

Herschel Reference Survey 

9.21 

8.62 

2.42 

8.39 

Herschel Reference Survey 

9.21 

8.57 

3.66 

8.47 

Herschel Reference Survey 

9.22 

8.69 

2.86 

8.78 

Herschel Reference Survey 

10.4 

9.05 

15.9 

9.66 

Herschel Reference Survey 

9.47 

8.68 

2.19 

8.43 

Herschel Reference Survey 

9.60 

8.84 

3.83 

8.33 

Herschel Reference Survey 

9.45 

8.84 

1.55 

<8.31 

Herschel Reference Survey 

9.17 

8.74 

7.54 

<8.16 

Herschel Reference Survey 

8.94 

8.92 

2.04 

<8.79 

Herschel Reference Survey 

9.27 

8.59 

3.63 

<8.56 

Herschel Reference Survey 

9.20 

8.82 

0.701 

<8.47 

Herschel Reference Survey 

9.23 

8.73 

1.61 

<7.98 

Herschel Reference Survey 

9.14 

8.60 

0.983 

<8.71 

Herschel Reference Survey 

11.48 

8.99 

4.29 

9.36 

LVL 

9.736 

8.36 

0.37 

7.85 

LVL 

10.86 

8.90 

2.45 

9.22 

LVL 

10.83 

9.31 

5.20 

9.12 

LVL 

9.485 

8.64 

0.13 

7.22 

LVL 

8.407 

8.01 

0.05 

8.14 

LVL 

8.450 

8.01 

0.07 

8.19 

LVL 

9.976 

8.76 

0.37 

8.35 

LVL 

10.74 

9.09 

0.97 

8.67 

LVL 

11.12 

9.04 

0.65 

8.66 

LVL 

11.10 

9.13 

3.10 

9.51 

LVL 

8.582 

8.23 

0.04 

7.76 

LVL 

9.722 

8.65 

0.17 

6.98 

LVL 

8.769 

8.25 

0.16 

7.04 

LVL 

9.296 

8.32 

0.31 

8.78 

LVL 

9.611 

8.80 

0.11 

7.21 

LVL 

11.01 

8.89 

1.91 

8.77 

LVL 

9.382 

8.31 

0.74 

7.92 

LVL 

9.376 

8.64 

0.28 

7.76 

LVL 

10.64 

9.01 

0.57 

8.59 

LVL 

10.45 

9.01 

0.47 

8.60 

LVL 

10.47 

9.05 

245 

10.75 

MS 

10.43 

8.95 

30 

10.44 

MS 

9.780 

8.73 

33 

10.85 

MS 

10.75 

8.73 

97 

10.70 

MS 


Table 2. (cont) Physical parameters (stellar mass, metallicity, SFR, and M(H 2 ) derived using the Wolfire et al. (2010) CO/H 2 conversion 
factor) for the galaxies in this work. 
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log M* 
[Mq] 

12+log(0/H) 

SFR 

[Mo/yr] 

M(H2) 

[M©] 

Source 

10.75 

8.84 

127 

11.07 

MS 

11.27 

8.61 

141 

11.07 

MS 

11.38 

8.96 

117 

11.02 

MS 

11.17 

8.98 

92 

10.80 

MS 

11.23 

9.01 

141 

11.17 

MS 

10.79 

8.75 

540 

11.54 

SMG 

11.20 

9.09 

810 

11.36 

SMG 

10.30 

8.75 

1070 

11.02 

SMG 

10.86 

8.74 

390 

11.48 

SMG 

10.59 

8.91 

680 

10.72 

SMG 

10.99 

9.07 

2890 

11.49 

SMG 

11.01 

9.09 

1260 

11.39 

SMG 

11.13 

8.86 

950 

11.07 

SMG 


Table 2. (cont) Physical parameters (stellar mass, metallicity, SFR, and M(H 2 ) derived using the Wolfire et al. (2010) CO/H 2 conversion 
factor) for the galaxies in this work. 
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